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ABSTRACT 



We calculate the statistical entropy of a quantum field with an arbitrary spin propagating 
on the spherical symmetric black hole background by using the brick wall formalism at 
higher orders in the WKB approximation. For general spins, we find that the correction 
to the standard Bekenstein-Hawking entropy depends logarithmically on the area of the 
horizon. Furthermore, we apply this analysis to the Schwarzschild and Schwarzschild- 
AdS(dS) black holes and discuss our results. 
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1 Introduction 



Thermodynamics and statistical mechanics of black holes are the most exciting and rapidly 
developing areas of black hole physics. The analogy between thermodynamics and black 
holes was originally introduced by Bekenstein (TJ [2] by assigning the area Ah of the 
black hole to the entropy Sbh- This correspondence was placed on the formal footing by 
Hawking's discovery that black holes radiate thermally with the characteristic temperature 

(T H ) PS, 

Hk ( c 

TH= 2^\k B ~ 

where k is the surface gravity of the black hole and h,c,kB denote the Planck constant, the 
speed of light and the Boltzmann constant, respectively Upon using the above expression 
for the temperature, the Bekenstein-Hawking entropy is given by 

with the Planck length, £ p = ^jG^h/c^. 

There are several approaches to derive the black hole entropy. In the Euclidean ap- 
proach [3 El EJ El E] , using the analytic continuation to the Euclidean sector and imposing 
the Matsubara period (3 = T^ 1 , the free energy and hence entropy have been obtained for 
the regular Euclidean solution of the vacuum Einstein equations. In the context of the 
string theory, the statistical entropy for extremal [10J [TTJ [121 H31 EH] and near extremal 
[T5| [16] black holes is determined by explicit counting of the microstates associated with 
the D — branes. Apart from these, the calculation of the black hole entropy and the 
corresponding quantum corrections have been studied by using various methods like- spin 
networks [T7"[[T5]. entanglement between the degrees of freedom separated by Killing hori- 
zons [121 1201 [21] , the conformal anomaly methods [221 E31 [21] and the quantum tunneling 
approach [231 [2S]. Very recently, the Euclidean approach mentioned above, has been 
implemented to derive the logarithmic corrections to Bekenstein-Hawking entropy for ex- 
tremal as well as non-extremal black holes in the various dimensions by incorporating the 
zero modes [27J. 

There is another efficient method, proposed by 't Hooft [2H], of computing the black 
hole entropy. This formalism, commonly known as brick wall model, is the semi-classical 
approach wherein the gravitational field (metric) of the black hole is treated classically, 
while the remaining degrees of freedom, leaving outside the horizon, are handled quantum 
mechanically. The entropy of these quantum degrees of freedom, calculated via statistical 
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mechanics, is identified with the entropy of the black hole. An important ingredient in 
the brick wall model is the boundary condition on the probe fields near the horizon. Any 
quantum fields near the black hole have a crucial property. Namely, the density of states of 
the quantum mechanical Hamiltonian of the fields blows out in the vicinity of the horizon. 
This results in the divergence of the statistical entropy [29j EDJ EH [32]. In the context 
of brick wall approach, this divergence can be controlled by putting a static spherically 
symmetric mirror near the horizon, at which the fields are required to satisfy Dirichlet 
or Neumann boundary conditions [28]. In other words, the distance between the horizon 
and the mirror (brick wall) behaves as an ultra violate cut-off. The canonical (statistical) 
entropy obtained by using this prescription is finite and agrees with Bekenstein-Hawking 
formula ([2]). This procedure also holds for any spacetime endowed with the horizon [53] . 
The status of the brick wall model have been elevated to formal level in [31] . They found 
that the expectation values of energy-momentum tensors for quantum fields in the ground 
state in the brick wall model matches exactly with the difference between Hartle-Hawking 
and Boulware states. 

All the methods of computing the statistical entropy of the black hole, in spite of 
their differences in the underlying assumptions and methodology, correctly reproduce the 
same leading order result. However, the quantum corrections are generally different for 
different approaches. For instance, the coefficient of the logarithmic corrections to the 
entropy obtained from fluctuations around the stable canonical ensemble [35] is different 
from the one obtained using tunneling formalism [26] [36] . It is then natural to seek for 
the extension of the brick wall model beyond the leading order. Such an extension of the 
brick wall approach up to sixth order in the WKB approximation has been reported in 
Ref. [37]. They showed that for spherically symmetric black holes, in four dimensions as 
well as in six dimensions, the corrections to the brick wall entropy can be expressed as 
F(Ah) \n(Ah/£p) where Ah is area of the horizon and the form of F(Ah) depends on the 
specific details of the black hole. 

Although the brick wall approach is implemented to obtain the entropy of a quantum 
scalar field in various background geometries and in diverse dimensions [38], relatively 
less attention has been given for its generalization to other type of fields, e.g, fermions, 
photons or gravitons [39]. Any computation of the entropy for a quantum field with 
generic spin degrees of freedom, even at the leading order, seems to be difficult [ID]. It 
is therefore worthwhile to analyze the higher order WKB corrections to the canonical 
entropy for arbitrary spins. In this work, by using the formalism given in Ref. [37] we 
calculate, up to second order in the WKB approximation, the statistical entropy of a 
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massless quantum field with the generic spin (scalar,vector,spinor etc.) propagating in 
the 3 + 1 dimensional spherically symmetric black hole background. For the general spins, 
we show that the functional dependence of entropy on the horizon area is the same for the 
leading as well as second order terms. However, the contribution to entropy coming from 
the second order term is appreciable compare to the leading order. The total entropy is 
obtained by adding the contributions from the leading and second order results. 

The paper is organized in the following manner. In the next section, we briefly outline 
the method for extending the brick wall entropy to higher order in the WKB approxima- 
tion and generalize it for arbitrary spins. Section 3 is devoted to explicit calculations of 
the leading order and second order brick wall entropies. Applications of this approach to 
specific black holes are exhibited in section 4. We summarise our results in section 5. 



2 Brick Wall approach for arbitrary spin 

We consider a 3 + 1 dimensional spherically symmetric black hole metri 



ds 2 = -g(r)dt 2 + ^— + r 2 {d6 2 + sin 2 Odcf 2 ) (3) 

g(r) 

where the form of the metric coefficient g(r) depends upon the specific black hole. The 
event horizon is defined by the condition g(r) = 0. For a spherically symmetric black 
hole, the Hawking temperature is given by 

Hk hg'(r h ) 

Z7T 47T 

where k is the surface gravity and the prime denotes differentiation with respect to the 
radial coordinate. We consider the massless minimally coupled field <3> p with an arbitrary 
spin \p\ (\p\ = 0, 1/2, 1 ■ • •), propagating in the background gravitational field ([3]). The 
equations of motion for $ p are given by [4TI 142] 

P D$ P = (5) 

where p O is the generalized d'Alembertian operator 

„□ = -jd 2 - 2pr (l - r -^j d t + (r 2 g)^ d r [{r 2 g)^d r ] 

1 cos 9 1 

+ — -d e (sm 9d e ) + 2ip-^-d^ + —^d\ - (p 2 cot 2 6-p). (6) 
sin U sin 6 sin a v 



1 Throughout the paper we shall work with the signature (—.+,+,+), and c = G = fcs = 1. 



3 



For the static and spherically symmetric background, we can consider the ansatz 



%{t^6,<P) = e-* Et l h R p {r)S, 
Substituting this in Eq.flSJ), we get the radial equation 



(r 2 g)- p 



d 
dr 



2 „\p+idRp{ r ) 



(r 2 g) 



dr 



+ 



r 2 u 2 



9 



A 



R P (r) 



2ipr 



r dg 
2g dr 



(7) 



R p {r) = 0, (8) 



while the spin weighted spherical harmonics Sf r 



m \ ' 



satisfy [13 



1 cos 6* 1 

— -<9 e (sin 6d e ) + 2ip—^ r -d <i> + —rzdl - (p 2 cot 2 9 - p) 
sinfl sin 9 sm 9 



with 



(9) 
(10) 



X = £(£ + l)-p(p+l) ; £> \p\ 

and E, I and m are the energy, orbital angular momentum and azimuthal angular momen- 
tum associated with the given normal mode, respectively. We can recast (jSJ) into more 
appropriate form by transforming the dependent variable R p {r) as 

1 



R p (r) = R p {r) 



{r 2 g)"* 
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After performing some algebra, we get 



E^_g\ 
h 2 r 2 



Rv(r) + 



2ipE 

h 



1 rg' 
g 2g 2 



R p (r) 



-(p + 1) 



— (P + !) + -r^(P ~ 1) + 7T + — 
rg Ag z 2g r z 



R p (r) 



(12) 



Our task is to solve Eq. (JT2"|) for each component^) of the spin \p\ and determine the 
number of radial modes of T P (E) of the field $ p with the energy less than E. The total 
number of modes T(E) is obtained by summing over all possible values of p. When we 
perform the sum, the terms proportional to ipE in Eq. (jl2p vanish. Therefore, we can 
write Eq. ffT2"l) as 



K(r) + 



Vlir) 

h 2 



V 2 {r) 



R{r) = 



(13) 



2 For a given value of spin \p\, there are 2p+l equivalent spin states, e.g for fermions we have \p\ = 1/2 
andp= 1/2,-1/2.) 
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with the definitions: 

Vi(r) 
V 2 (r) 



1 

9 

(P 



E 2 
-1) 



2\ 1/2 



— (P+ 1) + -r-z(p- 1) + — + 
rg Ag 2 2g r 2 



(14) 



In almost all cases (except for 1 + 1 dimensions) it is impossible to get an exact analytic 
solution of Eq. (TT5T) and the WKB approximation (at the leading order) is frequently used 
in the literature. We here compute the higher order WKB correction for a quantum field 
with a general spin \p\. 

We start by considering the WKB ansatz for R p {r) in Eq. lfiU]) . 

1 



R p (r) 



V 7 ^ 



-en 



lJdr'Q{r') 



(15) 



where 0(0 is an unknown phase. In order to analyze the higher order WKB approxima- 
tion, we expand 0(0 in powers of h 2 [37J HI] as 



Q{r) = Y,Q2,{r)h 2 \ 



(16) 



i=0 



After substituting Eqs. (113]) and (ITB]) into Eq. (TT5]) and collecting the equal powers of h, we 
get 

i2\ V2 



Oo(0 
Q 2 (r) 
Q 4 (r) 



±Vi(r) = ±- ( -E 2 

s V 

3Qg(0 Qg(r) 
8Ql(r) AQ 2 (r) 2Q (r 
501 (r) 1 



V 2 (r) 



+ 



20o(0 4Qg( r 

3Q 2 (0vy(0-0 2 (0vi"(Q 

41/ 3 (0 



(4o 2 (o^(o+o^(0) 



(17) 



for h 2 and ft 4 orders respectively. In fact, all the higher order functions can be expressed 
in terms of potential VI (0 an d V^(r). Upon using the semi-classical quantization scheme 
[28] together with the series expansion of Q{r), we write the expression for the density of 
states T(E) with energy less that E as 

oo 

T(E) = J2^(E). (18) 
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Here V2i(E) denotes the number of states at i th order and it is given by 




) 



2|p|+l 



L 




r 2i (£) 



E 



dr 



d£(2£+l)Q 2i (r,£,p). 



(19) 



P =i 



We have regularized our computations by introducing the ultra violate and infrared 
cut-off [28] such that the quantum field $ p satisfy: $ p (r < Vh + e) = $p(r > L) = with 
f C 1 and L ^> r^. There are few points which we would like to emphasize at this stage. 
Normal modes of the quantum field in the region close to horizon undergoes an infinite 
blue shift and eventually gives divergent contribution to the density of states. In fact, this 
feature is always present when we study the quantum field theory in any spacetime with 
horizon H5J H7J US]. The ultra violate cut-off introduced above precisely regularize 
this divergence. On the other hand, the density of states for the massless field at the 
large distance (asymptotic) gives divergent contribution to free energy. This divergence 
is cured by imposing the infrared cut-off. For the massive case the asymptotic density of 
states behave smoothly. 

Another point which we would like to mention is that, in Eq.( ll9p we have replaced the 
sum over I values by the corresponding integral. The maximum of the I integral is chosen 
such that Q (r,£,p) is real. Such an approximation of replacing sum by integration is 
valid for the leading as well as higher order WKB analysis. For the general spin, situation 
is slightly different. In that case, there exist a lower limit on the £ integration and it is 
related to the consistency of the eigenspectra of the spin weighted spherical harmonics 
JHJH9] (f° r more details, see [50] ). This modification complicates the higher order of the 
WKB analysis. 

Using the density of states (JT9~|) . we can compute the free energy for bosons(— ) and 
fermions (+) from their standard expressions, 



where (3 is the inverse of the Hawking temperature. The statistical entropy associated 
with the free energy is given by 




(20) 




(21) 
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3 Leading and Second order computations 



In this section, we shall compute the statistical entropy of the quantum field with spin 
\p\ propagating on the spherically symmetric black hole background at the leading and 
second order in the WKB approximation. The total entropy is obtain by combining the 
leading and second order results for the corresponding species. 
Leading order: 

From Eq. (fT9l) . we have the following leading order expression, 

g(r)Xh 2 ^ 1/2 



r (£) 



fm 



E 



dr / dt(2t+l) 



9{r) 



E 1 



(22) 



where A is given by Eq.ffTOl). It is convenient to work with variables A(r), E and G(A, E), 
which are defined as 



r 
E 2 



A 
E 

G(\,E) 

Using these, we rewrite the density of states as 

-L 



p(p+l)), 



~ ~\l/2 

E — X 



r (E) 



dr 



r 



9 2 (r) 



E 



dXG{X,E) 



(23) 
(24) 
(25) 

(26) 



The maximum of A is located where the function G(X,E) vanishes and the lower limit is 
dictated by the spin \p\. Using Eqs. liTUj) and (|25|) . we find 

g(r)H 2 



A, 



E ; A 



[\p\(\p\ + l)-p{p+l)} 



After performing the A integral and summing over p, we obtain 



r (£) 



(2\p\ + l)E 3 / 2 



j. — 

9 2 (r) 



Substituting this in Eq.( f20l and integrating over E, yields 



2 

45 



vr 3 (2b| 



14 /tt 3 (2|p| + 1) 



360 



h 3 (3 4 



r h +e 9 2 (r) 



r 

9 2 (r) 



(27) 
(28) 

(29) 
(30) 
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for bosons and fermions, respectively. 



Now, the radial integrations are evaluated by expanding the metric near the horizon, 



g(r) « g'(r h )(r - r h ) + 



9"(r h ) 



{r-r h f + 0{{r-r h ) 3 ). 



At the quadratic order in the metric expansion, we get 



7r 3 (2H + 1; 



1 

45 

7 1 V vr 3 (2H- 

8 45 V h 3 P A 



e 2 K 3 



+ 



rlg"(r h ) r h 



K 



Pk 3 



Ak 3 

r h9"(r h ) r h 



A In ^ 



2\ 1 



(31) 

(32) 
(33) 



An 3 k 2 

where we have used the coordinate invariant cut-off (proper distance of the brick wall 
from the horizon) 



(34) 



and a is the surface gravity given by Eq.(j4]). The canonical entropy computed from 
Eq.fFSTT) becomes 



v rlg"(r h ) r h K 



St 



rip) 



?2 



90e 
7 r 2 



rip) 



8 90e 



;-2 



f + ip)l 



In 

360 90 
r 2 h g"jr h ) _ run 
360 90 



(35) 
(36) 



Here f ± {p) = (2\p\ + l) is the spin degeneracy factor for fermions and bosons respectively. 
Our leading order results are in agreement with that of given in the literatures 



Second order analysis: 

We now calculate the density of states of the quantum field with energy less than E at 
the second order (i = 1) in the WKB approximation. Upon using Eqs. fl23|) . fl24|) and ( 125]) 
we can write Eq.( !T9|) as 



where 



I(E,r) 
J(E,r) 



T 2 (E) = — 



Th+e dr \gir) 



I{E,r) + J{E,r) 



An 



G(X,E) G 3 (\,E) 2 G 5 (A,E) 2 



E 

^ A, G(X,E) 



(37) 

(38) 
(39) 
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and we have denned 
Q?\r) 

g«(r) 



Q 



(2), 



9'{r) 
2r ' 

g'\r) 3g'(r) | 3g(r) | g"{r) 
8g(r) 4r 4r 2 8 
5 g' 2 {r) 5 g'( r ) 5 g(r) 



32 ^(r) 8 



Q 



(0), 



(p + 2) 



r 

<?'(r) 



0*(r) 



(p+1 



. g"( r ) 



2r 2 



4 



(40) 
(41) 

(42) 

(43) 



2r ~ &g{r) 

It is worthwhile to point out that Eqs. (j40j) .(141j ) and (I42p have the similar structures to 
those of given in Ref . [37] . While the additional piece Q^ir) depends explicitly on the 
spin orientation. This term would vanish for p = (scalar). Therefore, we expect that 
the A integration over the terms Q 2 > Q2 anc ^ Q2 should yield the similar forms given 
in Ref. [37J. However, the corresponding integral over gives completely different 
structure. We shall evaluate these integrals one by one. 

The definitions (124)) and (1251) allow us to write inverse powers of G(X, E) appearing 
in Eq. fl37j) as the partial differentials with respect to E. By successive applications of the 
Leibniz rule of integral calculus and by interchanging the order of A integration with E 
differentiation, we obtain 



I(E,r) 



(2\p\ + l)E l l 2 
AE 

+ 



2Q 



(o), 



±(° (- 



4Q 



E 2 



(i) 



16 ( 2 ) 



G(X, E) 3 ydE \2G(X,E) 
In getting the above expression, we have used 

2|p|+l 

^(E-A ) 3/2 (2^ + 3Ao) : 



E 2 
4G(A, E) 



(44) 



J X=E 



2E^ 2 (2\p\ + l] 



P =i 



2W+1 

(E-X f 2 (8E 2 + 12EX + 15X 2 ) = 8E 7 / 2 (2\p\ + 1) 



where A is given by Eq. fT27j) . 



Next, we substitute Q 2 i n Eq. (|39l and evaluate the integral to get 



2|p|+l 

J(E,r) = -2 Y,(E-X )V 2 p 
p=i 



(P + 2) 



g'(r) g' 2 (r) 

r p 

2r 8g(r) 



9{r) , g (r) 



(45) 
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As mentioned before, this term vanishes for \p\ = 0. For \p\ > 0, we can execute the sum 
by noting that 

2bl+i 

^P(E-X ) 1/2 = E^(l + \p\)(l + 2\p\), 
P =i 

2|p|+l 

Y,P 2 {E-\,fl 2 = -E 1 / 2 (l + |p|)(l + 2H)(3 + 4b|). 
p=i 6 

Using the above formulae and introducing the theta function, we write 

J(E,r) = -E 1/2 J(r)Q(\p\) (46) 

where 



J(r) = (l + b|)(l + 2|p|) 



(9 + 4H)^ + (3 + 4H)^i + (3 + 2H) 2 ^ r) ■ 9 " {r) 



3r v 1 u 12g(r) v 1 u 3r 2 2 

(47) 

Note that 0(|f>|) for \p\ = (scalar) vanishes, while it is unity for |p| ^ 0. 

We should emphasize that when A = X max — E (turning point) the function G(E, A) 
vanishes and all the derivatives with respect to E diverge. Consequently, integration of the 
terms like G _3 A and G~ 5 A 2 in Eq.( l38j) leads to the finite as well as diverging contributions. 
The divergent part is given by the last term in Eq. ( 144 p . This divergence is due to the fact 
that the WKB approximation is not valid as we move closer to the turning point. On the 
other hand, we observe that the result for J(E, r) contains only finite term. Thus, the 
structure of the divergent terms appearing for scalar particle [37] remain unchanged even 
for the generic spin. 

Then, the density of states V2(E) can be obtained by collecting the finite part of 
Eqs.( l4"4"|) and (j4"5"|) and substituting them into Eq. (157|) . 

E(2\p\ + 1) f L ^ r 2 |Y o _ (0 ), , Arii i), , , 16_ (2)( 



r a (J5) = V l h ' J dr— | \2QW(r) - 4Q^(r) + ^(r) ) 



-e(b|)(i + b|) 



(9 + A\p\)g'(r) , (3 + 4\p\) g' 2 (r) , 2(3 + 2\p\)g(r) , g"(r) 



3r 12 g(r) 3r 2 

Substituting the above expression in Eq.( l20l and integrating over E, yields 
(2|p| + l)i f L j_ r 2 /^ O (o), , , n( i), , « n(2) , 



F >- - --"w^ L, ir w) \ l 2er(r) - m + s « w J (49) 
_ e( |p|)(i + ipd (<i±iMMi) + e±ggp + w±m>m + m)y 
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(2\p\ + l)vr 



dv 

r h +e 9(r) 



2Q<?\r)-4Q£\r)+™Q<? ) (r) 



12/3 2 /i 

6(|p|)(l + \ p \) ( " (9 + 4|p|Mr) + (3 + 4|p|)g a (r) + 2(3 + 2\p\)g(r) + 



(50) 



3r 12#(r) 3r 2 ' 2 

Finally, by using the metric expansion ( 13 ip and performing the radial integration, we get 
(2|p| + 1)tt 



4r 2 



+e(b|)(i + b|) 



3eV(r A ) 
(3 + 4|p|) 
6 



2r fe g"(r h )r 2 h 
3 6#'(r fe ) 



(51) 



r 2 



xln^f 

(2b[+l>r 
12/3 2 n 

+e(b|)(i + b|) 



4r 2 



3eV(r h ) 



c2 9'(r h ) 



2r h g"{r h )rl 



In 



(52) 



(3 + 4|p|) / r 



6 



e 2 g'(r h ) 



j(7 + 4b|) 



2Ag'{r h ) 



(9 + 4|p|) 



x In 



Thus, the canonical entropy at the second order in the WKB approximation for the 
quantum field of spin |p| is given by 



So 



f-(p) 



r\ ( Kr h g"{r h )rl 



9e 2 



9 



+ 



72 



In 



(53) 



-e(b|)(i + b|) 



(3 + 4|p|) (r\ 



72 



^ - ^(7 + 4b|) + 



21 



288 



(9 + 4|p|) 



x In 



18e 



;2 



18 



g"(r h )rl 



144 



+e(b|)(i + b|) 



(3 + 4|p|) /r 2 



144 



x In 



In 



(54) 



-|^(7 + 4|p|) + 



ng"(r h ) 

576 



(9 + 4|p|) 



Eqs.( )53|) and ( 1541) constitute the main result. For the scalar field, our result is consistent 
with that given in Ref. [31]. The canonical entropies (I53p and (154]) contain the quadratic 
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as well as logarithmic divergent terms in the ultra violate regime of e — > 0. However, these 
divergences are even present at the leading order of the WKB approximations ( 135]) and 
( 136]) . Various approaches have been suggested in the literature for regularizing these diver- 
gences [2§] ED E2], and it was shown that they can be absorbed into the renormalization 
of the coupling constants appearing in the one-loop effective action [HU |52j [53] EH 155] . 

We now write the expression for the total canonical entropy up to the second order in 
the WKB approximation by adding the contributions from the leading (I35|36p and second 
( l53|54j) order expressions, 



s~- 



D 



Si 



(55) 



It is worthwhile to point out that the invariant cut-off e can be adjusted to match the 
leading order term with the standard Bekenstein-Hawking entropy (0). For the scalar 
field, the degeneracy factor f~(p = 0) is unity. Thus, by setting e 2 sc = (ll£p/907r) the 
standard expression for Sbh , at leading order, is obtained [37]. In general, the value 
of the invariant cut-off depends upon the type of the field. The massless spin- 1/2 (Weyl 



2 and in this case we set e^j 



(169£/1440tt). For 



fermions) we have f + (p = 1/2) 
the electromagnetic field (p = 1) we have only two polarizations corresponding to two 
helicities while the remaining one is unphysical gauge degree of freedom. Hence, the 
effective spin degeneracy factor for this case is f~(p = 1) — 1 = 2. The appropriate 
invariant cut-off length for electromagnetic field (photon) is ei h = (57^/907r). Thus, we 
write the corresponding expressions for the canonical entropy, 



s: 



s. 



ph 





K 


+ 9 - 


uf 


I0r h 




A h 


"471 




uf 


360 




A h 


"122 




Uf 


.Uf 





9 "{r h ) 
60 



An 
Air 

657 



In 



4320 
77 



g"(r h ) 



90 



9"{r h ) 



Ah_ 

167T 



hi 



Ah 

I 2 



167T 



.4, 



(56) 
(57) 
(58) 



for the massless scalar, the fermion and the electromagnetic field, respectively. 

The equations (153]) and (154]) represent the second order WKB correction to the lead- 
ing (zeroth) order WKB approximation (see Eqs (I35II36I) ) for bosons and fermions, re- 
spectively. It is possible to compute the higher order WKB corrections to the canonical 
entropies. Indeed, for scalar field, these corrections were presented in [27J. At the fourth 
order (i = 2 in eq. (119j) ). the density of states is T^(E) oc E^ 1 and therefore free energy 
becomes independent of f3 (i.e, scale invariant). Hence, the A th order contribution to the 
canonical entropy vanishes. For the sixth order (i = 3), density of states is T e (E) oc E~ z 
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and the corresponding expression for canonical entropy is given by 



S 6 _ 13892tt 2 r 2 
7~ " 45045 ¥ 



9vr 2 „, N 2 30vr 2 



g"( r h)r h + — Kr h 



77 u v " 77 



In ( J[ ) (59) 



where e; n in the above expression is due to the fact that density of states Tq(E) diverges 
in the infrared regime E — > [37J. Also, note that the UV divergent structure arising 
from 6 th order has the similar form as compare to the corresponding leading and second 
order terms. From the above analysis, it is natural to expect that the brick wall entropy 
does indeed provide corrections to the Bekenstein-Hawking entropy at all orders in WKB 
approximation. We hope that the generic structures for fourth and sixth order WKB 
corrections to the canonical entropy for the scalar field continue to be valid for the field 
with arbitrary spins. There are some important points which we would like to emphasize. 
As mentioned earlier, the WKB approximation arises when we try to solve the radial 
differential equation (I12p for the modes of the quantum field under consideration. In the 
standard brick wall formalism, existed in the literature, the canonical expressions for free 
energy and entropies were computed at the leading order in the WKB approximation. In 
our work, we obtain the possible modifications to the standard (leading order) result by 
taking into account the higher order WKB terms. This is different from considering the 
higher order loop corrections to the gravitational partition function. For instance, at the 
tree level the gravitational partition function reproduces standard Bekenstein-Hawking 
entropy [5]. While, at one — loop level, Bekenstein-Hawking entropy receives logarithmic 
corrections [2TJ 123 ETJ- The situation is slightly different in the brick wall formalism. 
Here, even at the leading order in WKB approximation, one can get the logarithmic 
correction [IB] (see also Eqs. ( 135| 136]) ). These logarithmic corrections appear when we 
expand the metric near the horizon upto the quadratic order in (r — 77J (see eq. (l3ip ) and 
this continues to be hold true irrespective of the orders in WKB approximation. However, 
the coefficient infront of the logarithmic term is different for different orders of the WKB 
approximation. Moreover, as we shall show below, for the Schwarzschild black hole, the 
leading order WKB result given in Eq.( )35|) exactly matches (including the coefficients) 
with the previous result given in [2Tj . 



4 Specific examples 

We apply the previous analysis for the Schwarzschild and Schwarzschild-AdS(dS) black 
holes and write corresponding expressions for the brick wall entropy up to the second 
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order in the WKB approximation. 



Schwarzschild black hole: 

The metric coefficient g(r) and the surface gravity k for the Schwarzschild black hole of 
mass M are given by 

g(r) = 1-^, (60) 
r 

« = ^- (61) 

where = 2M. Now we first substitute these expressions in the leading order expression 
(]35p for the brick wall entropy for scalare field and obtain 



S, 



(0)sc 9 q~2 



r_h_ 

360 V r\) 90 \2r h 



ln(^) =-L4- + ^ln[^). (62) 



This expression is in complete agreement with the corresponding result in [21]. The 
expression for second order canonical entropy can be obtain from Eq. (|53p and it is given 
by 

Sr 2)sc = A-^(^)- (63) 



(2)sc gp lg V g 

We now use Eq (}56"]) and write the expression for total entropy S 1 " = S^ sc + SZ\ ac interms 
of the horizon area A/j, as 

S „ = ^ l ln ^ (64 ) 



4£2 60 V ^ / 

Similarly, the total entropy for fermions and the Maxwell fields is given by 
+ A h 23 (A h \ A h 131 (A h \ 

Schwarzschild- AdS(dS) black hole: 

In this case the metric is given by 

2M 5A <A . o ( 2M kA 



ds 1 = ^1 - — -r 2 j dtr + ^1 - — -^j dr^ + r W (66) 

here A > is the cosmological constant and 5 = —1(1) for AdS (dS) cases, respectively. 
For S-AdS, there is only one horizon while for S-dS black hole we have event horizon (at 
r = Th) as well as the cosmological horizon (at r = r c ). 
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First we consider the S-AdS case and compute the canonical entropies. Comparing 
f )66|) for 5 = — 1 with eq.([3]), the metric coefficient and the surface gravity are given by 



9{r) - 
where satisfies 



r 



1 + (r 2 + rr h + r 2 h ) 



H 



1 

2^ 



2M = r h 1 



(67) 



(68) 



Substituting the above equations into Eqs. (l56|) -( 158|) . we obtain 



Ai 2 



Ah 
Ai; 



1 

60 4 



S~ = 



4^ 



131 



+ 



A 

SOtt" 
61 



144 240 



In 

AA fe 

7T 



c 2 

p 



A i 



-£\ ; S+ = 



wf 



Ai 2 p 



23 471 / AA h 
96 360 V 32tt 



In I 6 



Note that, unlike the Schwarzschild case, here the coefficients of logarithmic terms depend 
on the area. 

Now consider the Schwarzschild-deSitter black hole. We first compute the entropies 
at the event horizon which is located at r^. The metric coefficient and the surface gravity 
for the event horizon are given by 



9{r) 
and rh satisfy 



1 - 



A 



(r 2 + rr h + r 2 ) 



H 



2r h 



2M = r H [\ 



A 



(70) 



(71) 



Substituting the above equations into Eqs.f l56|) -f l58|) . we obtain 



s: 



s 



ph\h 



A h 


' 1 


A A 

307T 


Aif 


60 


A h 


"131 


61 


Aif 


144 


240 



111 

AA h 

7T 



Ah 

i 2 
1 p 



In 



■ 9 + I 

A h 



Ah 
Ai 2 p 



23 
96 



471 / AA 



360 V 32tt 



In 



i 2 



(69) 



(72) 



Next, we calculate the canonical entropies near the cosmological horizon which is 
located at << r = r c = \J\- We first consider the brick wall near the cosmological 
horizon with the proper distance e c . Note this proper distance is in general different from 
the one used to construct the brick wall near the event horizon. 
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In this case the surface gravity is given by k c = —^r c . Using this and following the 
similar procedure given above we arrive at 

c— I A 289 MAA /A c \ 

= 4£f + T35lT67j ln UJ- (?3) 

This gives us the entropies due to scalar, Weyl fermions and photons at the cosmological 
horizon. The total entropy for the S-dS black hole is then sum of the the entropies of 
the corresponding species at the event horizon and cosmological horizon. At the leading 
order our result is in agreement with the one given in [4TJj . 

We would like to mention that for since the S-dS black hole is endowed with the two 
horizons with different surface gravities and consequently two different Hawking temper- 
atures the quantum field between these two horizon is in thermodynamically unstable 
equilibrium. For this reason the computation of entropy for S-dS black hole is not of- 
ten discussed in the literature (but see e.g [56]). The derivation presented here should 
be supplemented by performing thorough analysis of quantum field theory on the S-dS 
spacetime |27J EEJ EH] • For example, it was shown that for the ordinary black holes horizon 
the ground state in the brick wall approach corresponds to the thermally excited state 
above the Boulware vacuum [31]. In other words, expectation values of operators calcu- 
lated for this state are the differences between expectation values in the Hartle-Hawking 
and Boulware states. One has to check the above criterion for the S-dS black hole also. 



5 Summary 

In this work, we have computed the canonical entropy for a massless quantum field with 
arbitrary spin propagating in 3 + 1 dimensional spherically symmetric black hole back- 
ground up to the second order in the WKB approximation. The generic structure of 
the leading and second order expressions for the free energy as well as entropy remains 
the same. However, the second order term contributes significantly to the entropy than 
the leading order. The total entropy up to the second order was obtained by combining 
the leading and second order results. The total entropy obtained in this manner con- 
tains quadratic as well as logarithmic divergent parts. These divergences can be cured by 
renormalizing the gravitational coupling constant [3"T||4"8]. However, the coefficient for the 
renormalized gravitational constant is found to be different for different species of matter 
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fields [5l]. Consequently, the proper cut-off distance of the brick wall from the horizon 
depends upon the type of field. Thus, by choosing the invariant cut-off appropriately, 
we were able to express the total entropy as a combination of the standard Bekenstein- 
Hawking entropy (Sbh) an d the logarithmic correction. The logarithmic contribution to 
the black hole entropy have been discovered earlier in several different approaches such as 
[27] [T7] IT8| [22| [26| |5T] . However, the coefficient of the logarithmic term is generally found 
to be different for different methods. In our case, the prefactor of the logarithmic term 
depends on the type of the field. For the scalar field, our result matches with that of given 
in Ref. [37]. It is important to note that these ultra violate divergences appearing in our 
analysis are independent of the order of the WKB approximation, as can be easily verified 
by comparing the leading and the second order expressions. We have also encountered 
the divergence in the evaluation of A integrals. This divergence occurs near the turning 
point of the WKB potential. However, the structure of these additional divergent terms 
still prevails irrespective of the additional spin degrees of freedom of the quantum field. 
Finally, we have applied our analysis to the Schwarzschild and Schwarzschild-AdS (dS) 
black holes and obtained the expression of the statistical entropy for the massless scalar, 
fermion and electromagnetic field, respectively. For Schwarzschild black hole, the brick 
wall entropy at the leading order in the WKB approximation (which also contain the 
logarithmic term) matches with known result [2T] . 

In the present work, we have restricted our computations up to the second order in 
the WKB approximation. The WKB analysis for the scalar field presented in [37] shows 
that the brick wall entropy receives no contribution at 4 th order while, at the sixth order, 
it receive significant contribution. We hope that these findings also hold for the quantum 
field with an arbitrary spin. As mentioned earlier these higher order WKB corrections 
are generally different from the loop corrections to the gravitational partition function. In 
this sense the WKB corrections which we are getting have different origin. These WKB 
corrections are need to be addressed in a sophisticated manner and are subject to future 
investigations. Another important aspect that we would like to explore is the thermo- 
dynamical stability of the black holes. In the tunneling mechanism, it was shown that 
the inclusion of the quantum corrections makes the black holes stable via phase transi- 
tion Rcf.[60j. It will be worthwhile to study the phase transitions and thermodynamic 
stability of the black holes by using the brick wall approach at higher orders in the WKB 
approximation. We would like to address these issues in near future. 
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